Evolution of density perturbations in double exponential quintessence models 
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In this work we investigate the evolution of matter density perturbations for quintessence models 
with a self-interaction potential that is a combination of exponentials. One of the models is based 
on the Einstein theory of gravity, while the other is based on the Brans-Dicke scalar tensor theory. 
We constrain the parameter space of the models using the determinations for the growth rate of 
perturbations derived from data of the 2-degree Field Galaxy Redshift Survey. 
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I. INTRODUCTION 

In the past few years, it has become apparent that 
the energy budget of our universe is dominated by an 
unknown component called " dark energy" . The WM AP 
table of "best" cosmological parameters for instance, 
gives a 0.73 ± 0.04 abundance for it, and a value of its 
equation of state uj < —0.78. In order to reheve some 
problems of the popular AC DM model (like the fine tun- 
ing issue) , a dynamical A-term has been proposed as rep- 
resentative of the dark energy. It's more popular version 
is a slowly rolling scalar field named quintessence. Many 
alternative cosmological models have been proposed, and 
indeed it is a challenge the work of ruling out all the 
"incorrect ones" on observational grounds. For instance 
many different potentials for these self interacting scalar 
fields (quintessence) have been proposed. However, it is 
obvious the importance of the observational exploration 
of the proposed cosmological models, and in this paper 
we give a further step in this direction. 

A variety of quintessence self-interaction potentials 
have been studied. Among them, a single exponential 
is the simplest case. This last model has two possible 
late-time attractors in the presence of a barotropic fluid: 
a scaling regime where the scalar field mimics the dynam- 
ics of the background fluid with a constant ratio between 
both energy densities, or an attractor solution dominated 
by the scalar field. Some of these models has been stud- 
ied in references 0, S 0, [E| • Given that single exponential 
potentials can lead to one of the above scaling solutions, 
then it should follow that a combination of exponentials 
should allow for a scenario where the universe can evolve 
through a radiation-matter regime (attractor 1) and, at 
some recent epoch, evolve into the scalar field dominated 
regime (attractor 2). For this reason the combination of 
exponentials represents an interesting alternative. Mini- 
mally coupled models with double exponential potentials 
are studied in [1, @, meanwhile, Brans-Dicke (BD) mod- 



els of quintessence with this kind of potential have been 
studied, for instance, in Q. 

The aim of this short paper is to (observationally) 
check models of the universe with potentials that are 
combination of exponentials: V = Vi exp(— ac/)) -|- 
V2 exp(— /30) (both in minimally coupled (Einstein) and 
BD theories), by considering another aspect of structure 
formation: the galaxy motions and clustering, i. e., the 
evolution of density perturbations in the Universe. In the 
past few years, observations of the large scale structure 
of the Universe have improved greatly. The development 
of fiber-fed spectrographs that can simultaneously mea- 
sure spectra of hundreds of galaxies has provided large 
redshift surveys such as the 2-degree Field Galaxy Red- 
shift Survey (2dFGRS) and the Sloan Digital Sky Survey 
(SDSS). In particular, the Anglo- Autralian Telescope of 
the 2dFGRS has obtained the redshift of a quarter mil- 
lion galaxies. This collaboration has produced abundant 
data and technical papers about galaxy motions and 
clustering, and we will refer to some of this, in particular 
their velocity to density comparisons. 

The paper is organized as follows: in section II we out- 
line the main characteristics of the models, in section III 
the main aspects of velocity to density comparisons are 
exposed and the equation for the growth of perturba- 
tions is solved, in section IV the observational check is 
presented and interpreted, while in section V conclusions 
are drawn. 



II. THE MODELS 

In this section we supply the main equations character- 
izing the dynamics of the models with double exponential 
potential of the form: 



V = Vie- 



V2 
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where Vi, V2, a and /3 are free constant parameters. We 
study separately models with minimal coupling (basi- 
cally Einstein gravity) and BD models written in Ein- 
stein frame variables. We adopt also throughout the pa- 
per units with 8ttG — c— 1 
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A. Einstein's Theory 



B. Brans-Dicke gravity 



We analyze flat Friedmann-Robertson- Walker (FRW) 
solutions to Einstein's theory with two fluids: a back- 
ground fluid of ordinary matter and a self-interacting 
scalar fluid that accounts for the dark energy component. 

If we set the following relationship between the 
barotropic index of the background (7) and the free pa- 
rameters of the above equation ([T|): 



/9 



37 



then it can be show that 
air) = { 



1 



2-3e 



sinh[/i£;(T-}-ro)]}3^-° 



(2) 



(3) 



is an exact solution of the Einstein's field equations for 
the above potential. Q . Where 



IJ'E = (87 - a 



2^ /7(2-3e) 
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and To is a constant of integration and the following time 
coordinate r has been introduced instead of the cosmic 
time t: dt — a~dT. 

This solution is very interesting since, the aforemen- 
tioned relationship between a and (3 leads always to a 
transition from a matter-dominated phase of the cosmic 
evolution at high redshift, into a (late time) dark energy 
dominated phase. This is, precisely, the kind of feature 
observational data suggests the cosmic evolution should 
share. Other parameters of observational interest are the 
Hubble expansion parameter: 
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COth[^£;(T - To)] (5) 



the matter density parameter: 



^mir) = (1 - e) [cosh[/^£;(T - To)]] 
and the equation of state (EOS) parameter: 

uj^t) = - 



3(l-f7™(T))' 



(6) 



(7) 



This solution depends only on three free parameters 
(7,a,£), where e is the density fl^{z) of dark energy 
in the early stages of the evolution (high redshift z ^ 
1). Nevertheless we will led with a reduce parameters 
space(Q;,£) because of we fix 7 = 1, meaning cold dark 
matter dominance at present. Using CMB[^, nucleosyn- 
thesis [1| and galaxy formation^ observations, the pa- 
rameters space can be constrained to be: < a < 1 & 
< e < 0.0450.1 



^ Note that SN la is not useful to constrain the free parameter 



Now we study flat FRW exact solutions to BD the- 
ory with two fluids: a background of ordinary matter 
and a self-interacting BD scalar field fluid accounting 
for the dark energy in the universe. In this case we are 
faced with two relevant frames (the Jordan frame and the 
Einstein frame), in which BD theory can be formulated. 
There has been discussion on whether these two frames 
are equivalent p^. It is not our aim to participate in 
this controversy and for practical reasons (simplicity of 
mathematical handling) we chose the Einstein frame. 

In this frame, if one assumes the following relationship^ 
between the free parameters a and /3 in ([T]) and 7: 

37 4-37 



/3 = ^ + ^^, (8) 



where lu is the BD coupling parameter; then 
a{r) = {y 



1 

is an exact solution with: 



sinh[/iB(r + ro)]}°(f*- 



(9) 



a{l3 — a) 



2-3£ 
3^2(1 -e) 



MB = ^^^W:,^7T^ (10) 



and 



2x/uj + I - a 



2. UJ + 



(11) 



The time coordinate r and the cosmic time in the Ein- 
stein frame t are related by the following expression: 

dt = a~dr. The bar notation means we are working 
in the Einstein frame of BD theory. The other interest- 
ing cosmological parameters are the Hubble expansion 
parameter: 



Hir) 



2-3£ 



-a(r)-'^coth[/Xi3(r-ro)] (12) 



/ 3n2(l -e) 
the matter density parameter 

a„(r) = n'il ~ £){cosh[/iB(r + ro)]}-2 (13) 
and the equation of state parameter 



W0(r) = -1 + 



3(l-a^(r)) 



(14) 



a, see Indeed, it is well known the controversy about the 
degeneracy of supernovae observations \VX Il3l [13 . Ha . Hsl Il7| . 
^ Recall that this relationship leads to a very desirable pattern for 
the cosmic evolution. See similar discussion under subsection A. 
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This cosmological model also depends on three free pa- 
rameters (7, a, s). Like in the former case we fix the value 
of 7 = 1 and executing a similar analysis with CMB, 
nucleosynthesis and galaxy formation observations, the 
parameter space can be constrained: < a < 0.37 & 
< £ < 0.045@.3. 



III. PERTURBATION GROWTH 

We will perform here a preliminary study of the evolu- 
tion of the mass density contrast {S = dp/p) in the mass 
distribution, modelled as a pressureless fluid, in linear 
perturbation theory. The evolution of quintessence den- 
sity contrast in not considered in this paper for simplicity. 
We understand, however, that these perturbations might 
influence perturbations of background matter [To[|. This 
method is based on Newtonian mechanics, that is bet- 
ter suited to the study of the development of structure 
such as galaxies and clusters of galaxies. This compu- 
tation requires that we be able to isolate a region small 
enough for the Newtonian gravitational potential energy 
and the relative velocities within the region to be small 
(non relativistic) [Tl| . The equation of time evolution of 
mass density contrast is 



TABLE I: This table show the characteristic parameters for 
each model 

parameters Einstein's Theory BD gravity 



X 



sinh[^E(r + ro)] sinh[/iB(r + ro)] 

10-3a^ a(£-2a)+4 
2{3-a^) a(/3-a) 



2(3-a^) a(0-a) 



1 6w"'(l-£) 
2{3-a^) a^(/3-a)^ 



Equation (fTS]) has two linearly independent solutions, 
the growing mode Sm+ and the decreasing mode Sm-, 
which can be expressed in terms of hypergeometric func- 
tions of second type 2-^1- We get 



oc 



4 + 4-^' -4 



d I 1 „ 

- + -0.1 + -g,-X^] 
4 4 2 ^ 



and 



2HS„i - inGpniSm = 0, 



(15) 



where the dot means derivative with respect to the co- 
moving time. In this equation the relative contribution 
of dark energy to the energy budget enters into the ex- 
pansion rate H. We shall consider this equation in the 
matter dominated era, when the radiation contribution 
is really negligible. The linear theory relates the peculiar 
velocity field v and the density contrast by [ll| 



---g,l--g,-X 



where 



Q, — 

4 4-^' 4 



1)2 +4e. 

For r <C 1 and r -C 1 we can write 



(19) 



47r J \x-y\-^ 
where the growth index / is defined as 



dln5„ 
rflna 



(16) 



(17) 



To solve the equation (|15p for the evolution or perturba- 
tions, it is useful to rewrite it in terms of suitable vari- 
ables, allowing some simplification 



(PSrn 
dX^ 



dS 

X^{X\ + d)-eS„, = 0, (18) 



where the new variable X, the parameters c,d and e are 
characteristic of model as follow: 



^ Note again that SN la observations are not useful to constrain 
the parameter space 



(20) 



(21) 



For determining the growth index of the perturbations 
we use the growing mode dm+ (I20p and substitute into 
(flT)) . It is well known the biasing effect in galaxy forma- 
tion, i. e.; the relative perturbations in the galaxy field 
and the matter field, on a point-by-point basis, are not 
equal: 



-(x)=fe^(x), 
n p 



(22) 



where n refers to the galaxy number density and b is the 
bias parameter. The parameter /3 = ^ relates the growth 
rate / of the perturbations (and hence the velocity field of 
the galaxy motions) with the density bias b. In this sense 
astrophysicists speak on a velocity /density comparison. 
Indeed, a compelling agreement is seen to exist between 
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FIG. 1: Parameter space for the model based on Einstein 
gravity. The parameter a is now constrained to a rather nar- 
row region. 



the velocity and density fields, which offers one possible 
test for the gravitational instability picture for the origin 
of structure [22*1 . 

The 2dFGRS has measured the position and the red- 
shift of a quarter million galaxies and from the analysis of 
the correlation function, determined the redshift distor- 
tion parameter (3 with the bias parameter h quantifying 
the difference between the galaxies and the dark haloes 
distribution. Usin g th e estimated /3 and the method em- 
ployed by [l^ [20l. |21|| to determine the bias h, one may 
estimate / = 0.51 ±0.11 at the survey effective depth 
z = 0.15 . As we can see the growing mode depend of 
the free parameters (a, e) as the growth index / to. Now 
we use this fact to additionally constrain the above men- 
tioned parameter space. 



IV. OBSERVATIONAL CHECK 

Applying equation P7|) . in the flat FRW solutions to 
Einstein's theory, does not give further constrain on e, 
but it does on a, as Figure [T] show. 

Assuming a flat universe, the parameter e — ri(^(oo), 
is the amount of dark energy in the very early universe 
{z ^ oo). It is known that an appreciable amount 
of dark energy at that epoch would imply an expan- 
sion fast enough to prevent the formation of structure 
at z ~ 3, but this parameter had been already con- 
strained in Ref.[5], so it is not problematic that now 
it has not been additionally constrained. In this case, 
the velocity/density comparison allows to locate a in a 
rather narrow region, thus acting as a stringent selector 
of attractor-like solutions of the kind 0, H, 0, [1| . Two 
things are to be said. Firstly, Figure [T] shows the joint 
region (a,e) satisfying the velocity/density comparison, 
for a alone we might have a wider variation. Secondly, 



FIG. 2: Parameter space for BD gravity obtained with the 
growth index analysis. 
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FIG. 3: Final parameter space for BD gravity. The shaded 
region obeys all of the constrains considered here. 



though narrow the interval, still we can speak of a as a 
selector of a class of solutions. 

Figure [3] shows the parameter space that was obtained 
for the model based on Brans-Dicke gravity. Taking ac- 
count this result and the former obtained in Q we con- 
struct the final parameter space which obey all constrains 
considered so far. 

As we can see in Figure |3] the parameters {a,e) were 
already considerably constrained from the original obser- 
vations. 



V. CONCLUSIONS 

In the models presented here we found that the pre- 
liminary study of perturbation growth is a good tool to 
constraint the parameter space in quintessence models 
with a self-interaction potential that is a combination 



5 



of exponentials. This kind of potential often appears in 
fundamental theories like the string or supersymmetry. 
Besides, it can accommodate a pattern of cosmic evo- 
lution that is characterized by a late time, dark energy 
dominated attractor. 

Research on the origins and evolution of the large-scale 
of the universe is one of the hottest topics in cosmology. 
In this work, we have used the relation between the pe- 
culiar velocity field of the galaxies, the growth rate of 
perturbations and the density bias in galaxy formation 
to make another step in the observational check of two 
quintessence models, resulting in a further constrain on 
the parameter space of the models. We consider that 
above discussion adds another argument in favor of the 



use of exponential potentials in quintessence cosmology. 
We plan in the oncoming future proceed this works us- 
ing another cosmological probes, like CMB, for instance. 
It could be desirable to consider perturbations of the 
quintessence field also. 
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